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Quantum walks, whose dynamics is prescribed by alternating unitary coin and shift operators, 
possess topological phases akin to those of Floquet topological insulators, driven by a time-periodic 
held. While there is ample theoretical work on topological phases of quantum walks where the coin 
operators are spin rotations, in experiments a different coin, the Hadamard operator is often used 
instead. This was the case in a recent photonic quantum walk experiment, where protected edge 
states were observed between two bulks whose topological invariants, as calculated by the standard 
theory, were the same. This hints at a hidden topological invariant in the Hadamard quantum walk. 
We establish a relation between the Hadamard and the spin rotation operator, which allows us to 
apply the recently developed theory of topological phases of quantum walks to the one-dimensional 
Hadamard quantum walk. The topological invariants we derive account for the edge state observed 
in the experiment, we thus reveal the hidden topological invariant of the one-dimensional Hadamard 
quantum walk. 

PACS numbers: 71.23.An, 03.65.Vf, 05.30.Rt, 78.67.Pt 


I. INTRODUCTION 


Topological insulators have attracted much attention 
from various branches of physics, due to their j-miique 
surface states predicted by topological invariantiUj. Al¬ 
though many materials, such as HgTe, Bi 2 Se 3 , etc., have 
been identified as being topological insulators, the nec¬ 
essary requirements on the intrinsic parameters, e.g., the 
spin-orbit interactions and internal magnetic fields, are 
hard to meet. One proposed way to overcome this diffi¬ 
culty, is to use Floquet topological insulators, i.e., to em¬ 
ploy a periodic drive to bring a material from a topolog¬ 
ically trivial phase to the topologically nontrivial one&cl. 
By altering the drive sequence, not only the p^apieters, 
but even the relevant symmetries of the systenJjcl can be 
tuned, offering a versatile route to topological insulators. 

Since the Floquet topological insulator is defined by 
a unitary time evolution operator, the spectral proper¬ 
ties of its effective Hamiltonian are characterized by the 
quasienergy, which has a 2tt periodicity, in natural units 
where the time is measured in units of the drive period 
and h = 1. As a consequence, in the presence of chiral 
or particle-hole symmetries, surface states at tt as well 
as 0 quasienergy should be considered. The number of 
TT quasienergy states is a novel topologically protected 
quantity, and is a corresponding novel bulk topological 
number. Accordingly, 1-dimensional chiral (or particle- 
hole) symmetric Floquet topological insulators are char¬ 
acterized by two topological invariants, ZxZ (or Z 2 XZ 2 ). 
Thus, the Floquet topological insulator provides richer 
physics compared with the time-independent topological 
insulator. 

Floquet topological insulators have already been real¬ 
ized in various experiments: by fabricating cawpled heli¬ 
cal waveguides for laser pulses in fused silicall3, by irra¬ 


diating the surface |©f a topological insulator with circu¬ 
larly polarized lightO, and by “shaking”, i.e., periodically 
modulating an optical lattic^with trapped cold atopis 
to realize the Haldane modeltd and Hofstadter models. 
However, none of these experiments has yet identified the 
unique tt quasienergy states so far. 


There is a promising way to realize Floquet topology 
ical insulators using discrete-time quantum walksEflE£l 
(quantum walks for short). The dynamics of a quan¬ 
tum walk is implemented by combining two fundamental 
operators: coin and shift operators, which change the in¬ 
ternal degree of freedom and the position of a walker, 
respectively. Recently, such quantum walks have been 
experimentally realized-iiu various systems, such, as cold 
atoms IS, trapped ionaiJllS, optical fiber ^psoESl, bulk 
opticsEl r-fwl integrated photonic circuitsE3. It has been 
clarifie(£3’Ej that the discrete-time quantum walk is an 
ideal platform to construct Floquet topological insulator 
because of the high tunability of relevant symmetriesErLl, 
and the parameters which are essential to establish non¬ 
trivial topological phases. Motivated by this work, stud¬ 
ies of t 
started 


pological phase of quantum walks have been 
, and their connection with the jeatanglement 
in these walks has also been investigatedE)!. Remark¬ 
ably, edge states originating in the non-trivial topological 
phase at zero and tt quasienergy have been observed in 
an e^eriinent on a one-dimensional photonic quantum 
walkEj. However, the bulk topological invariants pre¬ 
dicting these edge states were not identified: edge states 
were observed at an interface between two regions with 
the same topological number. 


In the present work, we identify the hidden topological 
invariants of the Hadamard quantum walk realized in the 
experiment of Ref. |^. We generali^ approach used 
for chiral symmetric quantum walksE3c3, and find a pair 
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of integers, i.e., a Z x Z topological invariant. 

This paper is organized as follows. In Sec. we de¬ 
fine the one-dimensional discrete-time Hadamard quan¬ 
tum walk and show how it is related to more commonly 
investigated quantum walks. In Sec. Ill we define a gen¬ 
eralization of chiral symmetry for the Hadamard quan¬ 
tum walk, and give the formulas for the corresponding 
chiral symmetric timeframes. In Sec. tV we calculate 
the topological invariants of the simple and the split-step 
Hadamard walks, and illustrate the consequences of these 
invariants, the topologically protected bound states, by 
numerical examples. In Sec^M, we apply this formalism 
to the setup realized in Ref. |32| , and demonstrate that the 
bound states observed in the experiment are predicted 
by our bulk topological invariants. Finally, Sec. ^ is 
devoted to discussions and conclusion. 


II. HADAMARD QUANTUM WALKS 


same matrix elements as the Hadamard operator up to 
the position of the minus sign. Although this seems like 
a small difference, it can have far reaching consequences, 
as we will show below. 

A shift operator is complementary to the rotation op¬ 
erator in that it changes the position of the walker in a 
way that depends on the value of the internal degree of 
freedom. We will use two shift operators, defined as 

:= ^ (|x ± l)(x| 0 |±)(±| + |x)(x| 0 |t)(t|); (6) 

X 

We consider two types of quantum walks, constructed 
from the Hadamard coin and the shift operators above. 
The simple Hadamard quantum walk is dehned via its 
timestep operator as 

Ua{0{x)) = S-S+Cii{9{x)). (7) 

The split-step Hadamard walk has the timestep operator 


In this paper we consider one-dimensional quantum 
walks where the walker has two internal states, denoted 
by 1+) := and |—) := (0,1)^. The wave function 

of the walker reads 

<S)\s), (1) 

s—± 

where a: € Z is the discrete position and t G N is the 
discrete time. The time evolution is generated by the 
unitary timestep operator as 

\'f{t)) = U*\'f{0)), (2) 

where the timestep operator U is composed of a sequence 
of coin operators C and shift operators S, to be defined 
below. 

A coin operator C acts on the internal state of the 
walker while leaving the position x unaffected, 

CH[d{x)\ := ^ |x)(a;| 0 'H[e{x)]. (3) 

X 

We take as coin operator the generalized Hadamard op¬ 
erator, 

with a parameter 9 depending on the position x. This 
can be expressed using the Pauli matrices 



and the identity matrix uo -=^ 2 , as 

^.[9] = (73 = cos 0(73 -|- sin0(71. (5) 

Most of the previous theoretical work used as coin opera¬ 
tor the rotation = cos 0(7o — i sin 0(72, which has the 


Ub[9i{x),92{x)] = .5_C'H[02(a:)]^+CH[0i(a;)]. (8) 

For both the simple and the split-step quantum walk, 
the effective Hamiltonian H is a useful tool to understand 
their long-time dynamics. It is defined from the unitary 
timestep operator by 

U = (9) 

Stationary states of a quantum walk are eigenstates of 
the timestep operator U, 

Um = X, = ( 10 ) 

Here the quasienergy e, the eigenenergy of the effective 
Hamiltonian, has 27r periodicity: due to unitarity of U, 
Xe takes its values from the unit circle on the complex 
plane. 

III. CHIRAL SYMMETRY OF HADAMARD 
QUANTUM WALKS 

For a one-dimensional Hamiltonian to possess topolog¬ 
ical phases, it needs to have some symmetry that-.links 
positive and negative energy states to each otherlTEl. We 
suggest an extension of the concept of chiral symmetry, 
and show that both Hadamard walks possess it. This will 
later allow us to describe the bulk topology and protected 
edge states of the Hadamard walks. 

A. Chiral symmetry at nonzero energy 

As a starting point we introduce the chiral symmetry 
at nonzero energy for a Hamiltonian. Consider a system 
of free fermions, with grand canonical Hamiltonian 

nm n 


( 11 ) 
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where the matrix of the single-particle Hamiltonian H is 
a continuous function of some system parameters denoted 
by ^ € S. This includes all parameters that are subject 
to disorder. The requirement for chiral symmetry of the 
Hamiltonian reads 

TH{OT = -H{0, ( 12 ) 

with a unitary chiral symmetry operator T = pl = T"^, 
that acts in each unit cell independently, and is indepen¬ 
dent of the disorder realization 

A global, fixed onsite potential ()) G R, that is not sub¬ 
ject to disorder, 

(13) 

obviously breaks chiral symmetry, since gM. commutes 
with any T instead of anticommuting. However, all it 
does is simply displace the energy of all states by 4>. Thus, 
if H hosts topologically protected bound states, so will 
H', the only difference is that they will be at energy 
instead of energy 0. 

The same discussion applies to periodically driven sys¬ 
tems. The requirement of Eq. (hm, translated for the 
time-evolution operator using EpT^), but allowing for a 
constant shift of quasienergy, reads 

TUiOT = (14) 


walks where the coin operator is a rotation, and where 
r := ^ |a;)(a;| 0 (Ti. (16) 

X 

To reveal the chiral symmetry of the Hadamard walks, 
we start by rewriting Eq. i) as 

n[e{x)] = (17) 

with X = ~4' = 7r/2. Importantly, x s-iid ^re global, 
fixed parameters. Now as for the factor it com¬ 

mutes with all operators, and so it will only shift the 
quasienergy. 

To deal with the factor in Eq. ([T^), we first 

make a few observations. It can be broken up into two 
parts, as 

g-*X<^3 _ g-ix/2-o-3 _ g-ix/2'0'3 

Then we notice that commutes with both S+ 

and S-, since 

S'± = ^ ( lx ± l)(x| 0 ^ ^ 

X 

Finally, we point out the relation 

. e-A/2-T3 . ^ gix/2 <T3 ^ ^g-*x/2 a3)-l^ (19) 


Importantly, not only T, but also 0 € R is here assumed 
not to be subject to disorder, i.e., independent of the 
parameters If Eq. (|^ holds, th^perator (^) has 
chiral symmetry (in the usual senseEUO) , and it may host 
topologically protected edge states at e = 0 or tt. If this is 
the case, the original timestep operator U{^) will have the 
same topologically protected end states at quasienergy 
e = 4>, respectively, e = tt -|- 0. In the following we will 
not write out the arguments ^ representing the effects of 
disorder explicitly. 


B. Chiral symmetry of Hadamard quantum walks 


We follow the ^thod developed by some of the au¬ 
thors of this workEZI to describe chiral symmetry of quan¬ 
tum walks, adapted to deal with chiral symmetry at finite 
quasienergy. A quantum walk has chiral symmetry at fi¬ 
nite quasienergy if its timestep operator can be split into 
two parts, conjugated inverses of each other: 

{7 = • TE-^r, (15) 


with ()) £ R. We can directly confirm chiral symmetry 
of the above quantum walk by substituting U and U~^ 
in Eq. ( p^ ) into left and right hand sides in Eq. (14), 
respectively, and using T”^ = T. In order to find such 
a decomposition, the starting time of the period can be 
shifted, i.e., the walk can be described in a different time- 
frame. This process is detailed in Ref. for quantum 


which will be useful to show chiral symmetry. 

Using the results of the previous paragraph, we rewrite 
the timestep operators of the Hadamard walks, Eqs. (j^) 
and (H) in a timeframe where the chiral symmetry 
at nonzero quasienergy is explicit. For the simple 


Hadamard walk this reads 

Ua = c-^^Fa ■ TFX^T-, (20a) 

Fa = (20b) 

TE^^r = S+ ■ (20c) 

while for the split-step Hadamard walk we find 

Ub = e-^^*FB ■ TFfi^T- (21a) 


Fg = g-*ei(3^)/2-<T2 5-_ . g-ix<T3 . g-*e2(x)/2.a2^ ( 21 b) 

TFg^T = e-*«2(U/2-<T2 . e-ix<^3 5-_^ . ^-^ex^)/2■a2^ ^21c) 

with, in both cases, x = ~4’ = '^/2- We remark that 
chiral symmetry of the simple and split-step Hadamard 
walks is preserved even when the parameter 9{x) of the 
coin operator in Eq. (i depends on the position x in a 
disordered way. 


IV. TOPOLOGICAL PHASES OF HADAMARD 
QUANTUM WALKS 

Having established chiral symmetry for the Hadamard 
quantum walks, we now determine the bulk topological 
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invariants, controlling the number of edge states in these 
walkscZrtaO. We will follow the procedure developed 
in Ref. which expresses the topological invariants as 
winding numbers of parts of the operator F from Eq. (H) 
between eigen^aces of the chiral symmetry operator E. 
In Appendix an alternative procedure developed in 
Ref. ^ is presented. 

To briefly summarize, Ref. ^ states that 

I r d 

vq =- / dk—\n detE+_(fc); (22a) 

zni dk 

1 d 

= -— / dk—h\ detE_ ik), (22b) 

Stti 7_7r dk 


bulk, with 9(x) = 0, in a chiral timeframe, as defined by 
Eq. (po|). The operator at quasimomentum k reads 

(26) 


The parameter y = 7r/2 shows up in two roles here. Eirst, 
it works as a global phase, /c-independent factor. This 
cannot change the winding number. The second role is 
as a displacement of k by y. This again does not change 
the winding number which is calculated by integrating 
over the whole k space. Thereby, we are free to set y = 0, 
when we substitute Eq. (^ ) into Eqs. {M) and (|^). We 
obtain 


where E_|_(fc) is the part of F in the quasimomentum 

space representation that maps from the subspace of the 
Hilbert space where E = — 1 (i.e., the eigenspace of E 
belonging to eigenvalue —1) to the E = +1 subspace, 

while F _is the part of F that acts in the E = — 1 

subspace. 

To adapt the results of Ref. ^ to the Hadamard quan¬ 
tum walks, we need t o tak e two things into account. 
First, Eqs. (2^) and ( 21b ) give the matrix of E in a 
basis where the chiral symmetry operator is not diago¬ 
nal, that is E = (Ti. In such a basis, i.e., whenever 


E = 


0 1 
1 0 


m = 


a(k) b(k) 
c(k) d{k) 


(23) 


for the parts of E necessary for the topological invariants 
we have 

2E+_ = (1 1) = (a - 6 + c - d), (24a) 

2E__ = (1 -1) = {a-b-c + d). (24b) 

Second, the Hadamard walks have chiral symmetry at 
finite quasienergy (j). Instead of topological invariants vq 
and we thus have invariants and i^Tr+cp, which read 


27ri dk 


(25a) 


= — [ dk-^\n (a(fc) - b{k) - c{k) + d{k)) . 
2 tti 7_7r dk 

(25b) 


We obtained this by substituting Eqs. ( pd]) into Eqs. (|22D , 
and omitting a factor of 1/2, which does not change the 
winding number. 


A. Simple Hadamard walk 


using the shorthand 

9 9 r- 9 ±ttI2 

se+= cos - ± sm - = V 2 sm- 

2 2 2 


(27a) 

(27b) 


(28) 


The invariant at quasienergy ±7r/2 is the winding number 
of a loop on the complex plane, centered at with 
radius se=p. In the case in which the radius is larger 
than the distance of the center from the origin, the loop 
encircles the origin, and we have a winding number of -1-1. 
In the opposite case, the winding number is 0. Therefore, 
to calculate the values of the winding numbers, we need 
to consider 


= 2 sin0. 


(29) 


For the winding numbers, the above considerations give 


, ^ / (0,1) if 0 < 0 < TT, 

(!Z_,,/2,I/+^/2) - I -7r<0<O. 


(30) 


The topological invariants are not defined when 0 = 0 or 
0 = TT. In these cases, the time-evolution operator reads 
Ua{9 = 7r/2±7r/2) = ±ie-dk+-^/‘i)'y3^ with a quasienergy 
spectrum that has no gaps. 


Numerical examples 

We illustrate the above results on the topology of the 
simple Hadamard walk, by showing examples where edge 
states appear near a boundary between two bulks with 
different topological invariants. To this end, we define 
two sets of 0(x): 

0, : {0_ = -V4, 0+ = +7r/4}, (31) 

Op : {0_ = +3^/4, 0+ = +^/4}, (32) 


We first consider the simple Hadamard quantum walk, 
as defined in Eq. (^. We take a translation invariant 


where 


0_ := 9{x < —1), 0+ := 9{x > 0). 
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FIG. 1. (Color online) Examples for bound states as signatures of a topological phase boundary in single step Hadamard walks. 
Top: Contour maps of the probability distribution P{x,t) va x — t plane for the quantum walks of (a-1) UA{0a), Eq. (^l|), 
and (a-2) Ua{0i3), Eq. (^). The topological numbers {y--n/ 2 ,i'+TT/ 2 ) for the negative and positive x regions are shown in the 
figures. The position dependence of the probability distribution P{x, t) ai t = 100 is shown in (b-1) and (b-2) for the single-step 
Hadamard walk UAiOa) and Ua{0p), respectively. The insets in (b-1) and (b-2) show eigenvalues A = of the corresponding 
time-evolution operator UA{da,p) of the main figure. The eigenvalues corresponding to edge states are distinguished by red 
crosses from those of bulk states (green thick arcs). 


According to Eq. (^) , we expect that topologically pro¬ 
tected edge states at quasienergies er = ±7r/2 appear near 
a: ~ 0 for UA{0a), because the two bulk regions have dif¬ 
ferent topological numbers, while no edge state should 
appear for Ua{Sp) because of the same topological num¬ 
bers in both regions. 

We numerically simulate the time evolution |'I'(t)) = 
UW^{t = 0)) up to t = 100 and calculate the probability 
distribution 

P(x,t):= ^ |((a:|®(s|)|vl/(t))|^ (33) 

The initial state is set to 

|vk(t = 0)) :=(!+)+*|-))/y2®|0). (34) 

Figures ^ (a-1) and (a-2) show the contour maps of 
P{x,t) in the x — t plane up to the time step t = 60 
for C/^(0„), Eq. (lU), and Ua{9p), Eq. ( |3^ respectively. 
On the one hand, Fig. |l] (a-1) clearly shows that the 
high probability amplitudes stably remain near the ori¬ 
gin where the topological numbers change. On the other 
hand. Fig. (a-2) exhibits low probability amplitudes 
near x = 0 as expected. 

The presence/absence of edge states is further high¬ 
lighted in Figs. 1^ (b-1) and (b-2) showing snapshots of 


the probability distribution P{x, t) at t = 100 for UA{da) 
and Ua{9p), respectively. We also numerically compute 
the spectrum of the timestep operators of the single- 
step Hadamard walks UA{9a) and Usidp), and show the 
eigenvalues in the insets. Here we consider 

the finite position space from —L to L — 1 with L = 100 
and impose the periodic boundary conditions to —L and 
L — 1. Thereby, we have two boundaries at a: = 0 and 
—T, where 9{x) is varied. If one of these boundaries 
hosts an edge state at quasienergy e, so must the other 
boundary: this gives an extra double degeneracy of the 
edge states. Although this degeneracy is lifted, because 
wavefunctions of the edge states at the same quasienergy 
but opposite edges overlap due to their exponential tails, 
this is the correction that is exponentially small in the 
system size, in our case, below the numerical accuracy 
~ 10“^®. Consistent with our theoretical prediction, the 
eigenvalues corresponding to edge states (red crosses) ap¬ 
pear at e = ±7r/2 for UA{da), while they are not there 
for Ua{9p). This illustrates the validity of Eq. (^). 


B. Split-step Hadamard walk 

We next consider the split-step Hadamard quantum 
walk, as defined in Eq. (^, with two coin operators 
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FIG. 2. (Color online) The phase diagram of the split-step 
quantum walk Ub{ 01 , 62 ) defined in Eq. (|^. Gapped phases 
are indexed by topological numbers Eqs. The 

red solid and blue dashed lines indicate closing of quasienergy 
gaps around e = 0 and tt, respectively. The symbols in the 
phase diagram indicate the parameters of the numerical ex¬ 
amples we consider, i.e., upper and lower triangles for 0 -^, 65 
in Eqs. an d f|38| ), circles for 0ei in Eq. (|4o|), and rectangles 
for 0e2 in Eq. ([4^)7^ 


Cii[0i,2{x)] applied during one timestep. We emphasize 
that this cuiantum walk has been realized in an optical 
experiments. We take a translation invariant bulk, with 
Oj{x) = 9j, for both j = 1,2, in a chiral timeframe, as 
defined by Eq. (H)- The operator Fb reads 


Again, the parameter x does not affect the winding num¬ 
bers, then we can set y = 0. 

The calculation of the topological invariants follows 
the same lines as for the simple Hadamard walk. We 
substitute Eq. into Eqs. ( p^ and (|^). We obtain 


^ {s9i+se2+ - e^’^se^-sg,,-) ; 


(36a) 


1 d 

^ ^ y ('Sei-'Ses-h + , 


(36b) 


using the shorthands defined in Eq. (p^, and bearing in 
mind that the quasienergies are displaced by tt instead of 
7r/2 as for the simple Hadamard walk. 

Using the same logic as for the simple Hadamard walk 
for the winding numbers gives us 


Vtt = (sign[—sin0i — sin02] + 1) /2; (37a) 

vq = (sign[sin0i — sin02] + 1) /2. (37b) 


The phase diagram in Eig. p| of the topological num¬ 
bers vq and Vtt at quasienergies e = 0,7r, respectively, 
Oft/B(0l,02). 


Numerical examples 

We illustrate the topological properties of the split-step 
Hadamard walk using two parameter sets, 

: {01_ = 02- = -^/4; 01+ = 02+ = +7r/4}; (38) 

Gs : {01- = -02- = -7r/4; 0i+ = 02+ = +7r/4}, (39) 

where 

01(2)- •= 0 i (2)(3^ < ~1)) ^1(2)+ ■= 0 i (2)(3^ ^ 0). 

Because only 0 = ±7r/4 is employed, the quantum walk 
UB{G-y, 5 ) would be called the split-step Hadamard walk. 
As indicated by the triangles in Eig. ||, the parameters of 
the sets G^^s are located on the red solid and blue dashed 
lines indicating the quasienergy gap closing around £ = 0 
and TT, respectively. In case of G^, the quasienergy gap 
around e = 0 vanishes, while the other gap around e = tt 
is still open. Since the topological numbers differ in 
the positive and negative x regions for Ub{G^), the edge 
states should exist. However, in case of Gs, parameters 
of the negative (positive) x region locates on the blue 
dashed (red solid) line. This results in no more energy 
gaps in the whole system. Then, we predict no edge 
states for Ub{Gs)- 

We confirm these predictions from the phase diagram 
by numerical simulations as shown in Eig. |^. In case of 
C7_B(0-y), we confirm the edge state at the quasienergy 
e = TT as well as the gap closing around £ = 0. We 
also numerically confirm that no energy gaps emerge for 
Ub{Gs), and then no edge states. 


V. INTERPRETATIONS OF HIDDEN 
TOPOLOGICAL INVARIANTS IN EXPERIMENT 

Finally, we resolve the hidden topological invariant 
found in the photonic quantum walk experiment in Ref. 
[ 3 ^ where the time-evolution operator is the one given in 
Eq. (^). We consider two par^eter sets which are also 
investigated in the experimenfEj: 

0ei : (01- = 0, 02- = +7r/4; 

01+ = +77r/I6, 02+ = +^/4}, (40) 

Ge2 : (01- = -5^/16, 02- = 0; 

01+ = +I37r/I6, 02+ = 0}, (41) 

In Ref. 1^, it is reported that the parameter smces with 
the topological numbers (1, 0) and (0,1) in Fig. g have the 
topological number 0 and the regions with the topological 
number (0,0) and (1,1) in Fig. || have the topological 
number 1. Thereby, in the case of G^i, the topological 
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FIG. 3. (Color online) Examples of bound states as signatures of a topological phase transition in split-step Hadamard walks. 
Top: Contour maps of the probability distribution P(x,t) in the x — t plane for the quantum walk (a-1) UsiO-y), Eq. ( |38| ) and 
(a-2) Ub{0s), Eq. (^ol). The topological numbers {vo^vj) for the negative and positive x regions are shown in the figures. A 
letter “N” means that the topological number for the corresponding quasienergy cannot be defined because of the gap closing, 
(b) The position dependence of the probability distribution P{x, t) a,t t — 100. The eigenvalues of the timestep operators are 
shown in the insets as in Figs. (b-1) and (b2). 


numbers of the positive and negative x regions differ by 
1, and then edge states are expected. In the case of 9f.2, 
the topological numbers in both regions are zero, and 
then edge states are not expected. However, edge states 
are observed in both cases in the experiment. 

In Fig. H (a) and (b), we show the probability distribu¬ 
tion P{x, t) att = 100 of the walks UB{dei) and UB(de 2 ), 
respectively, starting from the initial state 

\ip{t = 0)) := |-h) (g) |0), 

which also coincides with the experimental. The corre¬ 
sponding eigenvalues of UB{dei) and UB{de 2 ) are shown 
in the insets. We confirm that the split-step walks 
Usidei) and UB{de 2 ) exhibit edge states at £ = 0 and 
£ = 0, TT, respectively. 

Now we look at the phase diagram in Fig. The 
phase diagram predicts that the quantum walk UB{9ei) 
should have an edge state at £ = 0 because the regions 
for a; < — 1 and a: > 0 have the topological numbers 
= (0,0) and (1,0), respectively. Furthermore, 
in the case of f7s(0e2), the edge states should appear at 
quasienergies £ = 0, tt because the regions for x < —1 and 
X > 0 have the topological numbers (0,1) and (1,0), re¬ 
spectively. These theoretical results are completely con¬ 
sistent with the numeri(^ results in Fig. and observa¬ 
tions in the experimentEHI. Thereby, we have succeeded 
in explaining the hidden topological invariant found in 


the experiment by the phase diagram Fig. || which is de¬ 
rived by establishing the relation between the rotation 
and Hadamard matrices. 

VI. DISCUSSIONS AND CONCLUSION 

Comparing Fig. (a-1) and (a-2) with Fig. (a), we 
notice that the the contour maps of the probability distri¬ 
bution P(x, t) in the former ones are sparser and checker¬ 
board like. The origin of this is a sublattice symmetry of 
the time-evolution operator defined as 

TsUTs = -U, (42) 

where 

Fs := ( |a;)(a;|- Y |a;)( 2 ;|) O cro- 

aGeven a:Godd 

The single step walk Ua retains sublattice symmetry, 
while the split-step walk Ub generally does not. This 
symmetry constrains the walker at every timestep to hop 
from the even sublattice (x even) to the odd sublattice 
(x odd), leading to the checkerboard pattern of P(x,t). 

If a quantum walk has sublattice symmetry, as defined 
by Eq. ( ^ , every eigenstate at quasienepjy £ must have a 
sublattice partner at quasienergy £-|-7rEj’E3. This explains 


































position X 

FIG. 4. (Color online) The position dependence of the prob¬ 
ability distribution P{x, t) at t — 100 is shown in (a) and (b) 
for the split-step walks UB{dei) and [/s( 0 e 2 ), respectively. 
The topological numbers {vo,Vn) for the negative and posi¬ 
tive X regions are shown in the hgures. The meaning of insets 
is the same as those in Figs. ^ (b-1) and (b2). 


where r] = x = —(j) = 'K/2 are global and fixed parame¬ 
ters. Thereby, cti also does not break chiral symmetry. 

In Ref. two topological invariants, and Q'^ are 
defined that are associated with a boundary between two 
bulks. These invariants are nothing but the topologically 
protected number of 0 and tt quasienergy edge states lo¬ 
calized at the boundary. The topological invariants vq 
and z/tt we define in this paper are associated with the 
bulks. By the arguments detailed in Ref. the change 
in vq as we cross from the left bulk to the right pre¬ 
dicts the value of (,Q^)- this is the bulk-boundary 
correspondence for one-dimensional Hadamard quantum 
walks. 

The formulas for the topological invariants in this pa¬ 
per assumed translation invariance in the bulk, so that 
the quasi-momentum fc is a good quantum number. For 
disordered quantum walks, e.g., where 0{x) is a ran¬ 
dom function of position, there are alternative formu¬ 
lations of the tnpnlngical invariants, based on the scat¬ 
tering matricesEStJ. These can be applied to disordered 
Hadamard walks using the mapping we presented in this 
paper. 

In summary, we have studied the topological phases of 
the one-dimensional Hadamard quantum walk. We have 
generalized the definition of chiral symmetry, and pro¬ 
vided a sufficient requirement for quantum walks to obey 
this symmetry, in Eq. ([T^. Employing the generalized 
definition, one-dimensional Hadamard quantum walks 
have chiral symmetry, and the corresponding topologi¬ 
cal invariants, which characterize the topological phases, 
can be calculated. We have used this result to reveal the 
topological invariants behind a recent photonic quantum 
walk experimenlEj. Our results add to the growing body 
of knowledge on the topological phases of quantum walks 
and Floquet topological insulators. 


why the single step Hadamard walk has edge states at e = 
± 7 r /2 appearing simultaneously. We note that, due to the 
specific value of 6*2 = 0 in the parameter set 0 e 2 , the split- 
step walk UB{de 2 ) has sublattice symmetry. Therefore, 
the edge states appearing at e = 0, tt in the inset of Fig. 
^ (b) are sublattice symmetric partners of each other. 

The results of this manuscript can be applied straight¬ 
forwardly to quantum walks where the generalized 
Hadamard coin is replaced by operators with 

j = 1, 2, instead of j = 3. In the case of j = 2, we again 
obtain a chiral symmetric walk: in analogy with Eq. (17), 
we have 




(43) 


where t] = —(p = tt/2. Thereby, CT 2 shifts the angle of the 
rotation matrix and the quasienergy by 7r/2, but does not 
affect chiral symmetry. The other case, j = 1, is easily 
understood because ai = Using Eqs. (0) and 

(^, we obtain 


g-jeo-2^^ _ ^-i4>g-i{S+r])a2 . g-ixo-3 


(44) 
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Appendix A: Alternative calcnlation of topological 
nnmbers 


In Sec. IV, we calculate the topological numbers of 
the quantum walks based on the method developed in 
Ref. Here, for completeness, we present the calcula¬ 
tion following the method of Ref. The two methods 
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are equivalent, we present this derivation for pedagogical 
reasons. 


1. Single-step Hadamard walk 


First, w e fo cus on the single-step Hadamard walk. 
From Eq. (20a), the chiral symmetric form of the single- 
step Hadamard walk is written as 


Ua = • Fa ■ FF^^F, (Al) 

_ g-i6»(x)/2-cr2 . g-jx/2-o'3^_ 

FF^^F = S+ ■ 

where x = —(/> = 7 r/ 2 . As explained in Ref. ^ if a 
time-evolution operator of the quantum walk U has chi¬ 
ral symmetry, another chiral symmetric time-evolution 
operator U' can be identified only at the different “time- 
frame”. The one for Ua^ thus is given by 


• FF^^F • Fa- (A2) 


Since the phase (j) only shifts the quasienergy, we proceed 
in our calculations by setting ^ = 0 on the above equa¬ 
tions during the calculation and shift the quasienergy by 
tt/2 at the end. In the momentum representation (by 
assuming the constant 0 ), we have 


UA{k) = e 


_ -ie/2-a2 


,-i(fc+x) 


0 


^- 16 / 2-02 


Q gik{k+x) 

= cos(fe + x)ceCTo - ickSe(J 2 - isin(fe + x)cr 3 , 


and 


U'A{k) 


/ ^-i(k+x) 

V 0 


0) „ ('i 


e 


0 \ 

^(k+x) J 


cos{k + x)cecro - - i sm{k + x)c0(T3- 


Here, we use the shorthands 


ce = cos( 0 ), se = sin( 0 ). 


The factor x only shifts the momentum k, which does 
not change the topological number. Then, we set x = 0 
in the following. 

Applying a unitary transform so that the time evo¬ 
lution operators have chiral symmetry in the basis that 
the chiral symmetry operator is diagonal, i.e., F = 173 , 
we have 

UA{k) = 

and 

= CkCgao + iskCgai — isga^- 

Since the coefficients of the (Tq term of UA{k) and lJ'A{k) 
are the same, both operators have a common eigenvalue 


Aa,± = 6='=“'^, sin(wA) = \/l - (cfcce)2 > 0. 


The corresponding eigenvectors \tp±) and \4’±) of UA{k) 
and U'^{k), respectively, also have the similar structures 




1 

71 V 1 y 




but with different phase factors 


1 

71 


(A3) 


F'^Aik) = (-CkSe + isk)/ sin(wA), (A4a) 

F'-f’A^k) ^ +iskCe)/sin(uJA)- (A4b) 


The winding number is defined through the Berry phase. 


1 


dk{'il)\d/dk\’4)). 


(A5) 


Substituting eigenvectors in Eq. ([^), the winding num¬ 
bers V and v' oiU a and C/^, respectively, become 


V = 


1 


dipA{k), 



Thereby, the winding numbers are determined from the 
trace of <fA{k) and ip'A(k) as k is changed from 0 to 27r. 
Considering Eq. ([^), this gives the following results, 

_ J -1 (0 < 6» < tt) 

^ \ 1 (-TT < 0 < 0) ’ 


and 


F = 0. 


Einally by using a formula derived in Ref. 27 in order to 
calculate the topological numbers for quasienergies e = (j) 
and TT + <h 


V + v V -V 

we obtain 

f (- 1 / 2 , + 1 / 2 ) for 0 < 6 » < TT 

(>2-^/2, >a+^/2) I (+1/2,-1/2) for- 7 r< 6 l< 0 . 

(A7) 

Since the global shift of topological numbers does not 
alter the argument of the bulk-edge correspondence, we 
confirm the consistent result with Eq. (jsOlin Sec. 0 by 
shifting numbers in the right side of Eq. (|A7D by 1/2. 


2. Split-step Hadamard walk 

In the case of the split-step Hadamard walk, we have 
the following two chiral symmetric time-evolution oper¬ 
ators 


Ub = • Fb ■ FF-iF, (A 8 ) 

U'B = e-^^^-rFB^r-FB, (A9) 
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where 


The eigenvalue oiUsik) is 


and 


Pg = ^-iSi{x)/ 2 -a 2 g_ . g-ixcrs . g-*62(a;)/2-(T2 ^ 

B ' ’ 

with X = ~4‘ = '^/2- Again we set </> = 0 of ?7 b and 
and shift the quasienergy by tt at the end of the 
calculation. We derive the time evolution operators in 
the momentum space representation as 

Usik) = [cos(fc + 2 x)ce 2 cei - se^sejcro 

- i[cos(A: + 2 x)ce 2 Sei + se 2 Cei]o -2 

- i sin(fc + 2x)c6i2cr3, (AlO) 

U'g{k) = [cos(/c + 2x)cg^ce^ - se^sgjcro 

- i[cos(fc + 2 x)cfliSe 2 + se^ce 2 ]o -2 

- ism(k + 2x)c6(ifT3. 

Comparing the above two equations, we notice that 
Usik) and U'g{k) are identical only by switching 9i 
and 02- This means that results for 17s(fc) are imme¬ 
diately applied to those for Ug{k) by switching 9i and 
02- Thereby, we present calculations only for Usik) here¬ 
after. 

Similarly to the si ngle-s tep Hadamard walk case, we 
can set x = 0 in Eq. ( AlO ) and apply the unitary trans¬ 
formation, we have 

UB{k) = 

[cOs(/c)c02C0i 
-I- i shi{k)ce^(Ji 
- i[cos{k)ce2S0^ + se,,ce^]a2. 


Ab.± = e±“^ 


sin(wB) = \/l - [cos(fc)c6(2Cei - > 0, 

and the corresponding eigenvector is 

I, , 1 (pie'‘VB(k)\ 

~ 72 V 1 / ’ 

^iVB{k) ^ -[cos(fc)ce 2 Sei + se^ce^]+ism.{k)cg^ 
sin(w_B) 


(All) 


Substituting Eq. (|All[) into Eq. (|A5[) , the winding 
number v oHJb is summarized as follows: when 


si7(0i) — si7(02) > 0, 


-1 (0 < di < tt) 

1 (-TT < 6»i < 0) 


(A12) 


(A13) 


otherwise 


v = Q. (A14) 

As we mentioned, the winding number v' of Ug{k) is 
given by switching 0i and 02 in the above results for v. 

Finally, we obtain the consistent phase diagram with 
that of Fig. H in Sec. 0 by substituting z/ and ly' oi Ub 
and Ug, respectively, into Eq. (A6), taking account of the 
quasienergy shift by tt, and 1/2 shift of the topological 
numbers. 
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